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Abstract 

We derive the wave equation obeyed by electromagnetic fields in curved space- 
time. We find that there are Riemann and Ricci curvature coupling terms to the 
photon polarisation which result in a polarisation dependent deviation of the pho- 
ton trajectories from null geodesies. Photons are found to have an effective mass 
in an external gravitational field and their velocity in an inertial frame is in general 
less than c. A physically relevant consequence of the analysis is that the curvature 
corrections to the propagation of electromagnetic radiation (in a homogenous and 
isotropic spacetime) keep the velocities subluminal provided the strong energy condi- 
tion is satisfied. We further show that the claims of superluminal velocities in higher 
derivative gravity theories are erroneous and arise due to the neglect of Riemann and 
Ricci coupling terms in the wave equation, of Einstein gravity. 



A standard result of Einstein's gravity is that the trajectories of all massless 
particles are null geodesies. A question worth examining is whether there is a devi- 
ation from the null geodesies if the particles have a spin i.e. due to the interaction 
of spin with the Riemann and Ricci curvatures of the gravitational fields. In this 
paper we have studied this question for the case of photon propagation in a curved 
background. Starting with the action y/^gF^ F^ v for electromagnetic fields in a 
gravitational field, we derive the wave equation for electromagnetic field tensor F^ u , 
which turns out to be of the form (Eddington [1] , Noonan [2]) 

W M F„ A + R^xF"* - R\F vp + R\F Xp = (1) 

We see that the photon propagation depends upon the coupling between the Riemann 
and Ricci curvatures and the photon polarisation. This leads to a deviation of the 
photon trajectories from the null geodesic by amounts proportional to the Riemann 
and Ricci curvatures. The photon trajectories in the geometrical optics limit are 
described by the following generalisation of the geodesic equation 

^ + T ^H7 = 2 V ° ^ ^ + - x fa (2) 

The nonzero right hand side of the modified geodesic equation (2) implies that there 
is a polarisation dependence in the gravitational red shift, bending, and Shapiro time 
delay of light even in classical Einsteins gravity. 

A curious phenomena discovered by Drummond and Hathrell [4] is that in 
higher derivative gravity which arises by QED radiative corrections, the photon ve- 
locity in local inertial frame can exceed the velocity of light in the Minkowski space. 
Due to the coupling of the electromagnetic fields with the Riemann and Ricci ten- 
sors in the Lagrangian, it was claimed that the photon velocity in the Schwarzschild, 
Robertson- Walker, gravitational wave and deSitter backgrounds is larger than the 
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flat space velocity c. This result was extended by Daniels and Shore to charged [5] 
and rotating blackholes [6] with the same conclusions. Latorre et al [7] have shown a 
universal relation between the velocity shift of photons to the energy density which 
generates the background metric, and Shore [8] has related the velocity shift to the 
coefficients of conformal anomaly. Lafrance and Myers [9] interpret these results as 
the breakdown of the Equivalence principle in higher derivative gravity and Dolgov 
and Khriptovich [10,11] derive this result from field theoretic dispersion relations. 
Finally Mende [12] has proposed this effect as a test for string theories of quantum 
gravity. 

In this paper we show that claims of superluminal photon velocity are due to 
the neglect of the Riemann coupling terms in the wave equation which arises from the 
minimal F^F^ Lagrangian. We find that in Einstein's gravity the photon velocity 
in a Schwarzschild blackhole metric and in the Friedman-Robertson- Walker metric is 
less than c. In other words the photon trajectories are always inside the null cone. 
We find that for the photon trajectories not to go out of the null cone the background 
matter should satisfy the strong energy condition p > 3p. This provides us with an 
answer to the question raised by Zeldovich and Novikov [13] - What law of physics 
would be violated if the strong energy condition is not satisfied ? Our answer is 
that special relativity in the free fall inertial- frame demands that the strong energy 
condition be satisfied. 

We also derive the wave in higher derivative gravity and show that Riemann 
coupling terms in the lagrangian of the higher derivative gravity are always smaller 
in magnitude than the Riemann term that already exists in Einstein's gravity. This 
analysis shows that the photon velocity does not exceeds c, even by the inclusion of the 
radiative correction higher derivative terms in the Lagrangian contrary to the claims 
[4-12]. The absence of superluminal propagation is not dependent on the eikonal 
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ansatz but is a consequence of the fact that the wave-equation in for the electromag- 
netic fields in Einstein and higher derivative gravity is hyperbolic. Consequently the 
field solution at a given point can only depend upon sources inside the past null- 
cone. In references [4-12] the vector potential solution is assumed to be of the eikonal 
form, and it is shown that the vector potential propagates outside the null-cone. The 
eikonal approximation for the vector potential does not always correctly describe the 
propagation of electromagnetic waves since the vector potential is not gauge invariant 
and it can be non-zero even in the acausal regions where the electromagnetic field is 
zero. 

The interaction of electromagnetic fields with gravity is given by the action 

S = J d'x^F^ (3) 

(We use the convention c = 1, signature —2 and Greek letters denote spacetime 
indices — 3). From (3) we obtain the equations of motion 

VF^ = (4) 

Equation (4) with the Bianchi identity 

V^F, A + V„F^ + V A F M „ = (5) 

gives the Maxwell's equation in curved background. Operating on equation (4) by 
the Va and using the Bianchi identity (5) the wave equation may be obtained as 

VV,^ + Va(V^) + [V^,V„]F v 

(6) 

- [V A ,V]F^ = 

The second term vanishes owing to (4). Using the identity for the commutator of 
covariant derivatives: 

[V", V„] F QjU = R pau ,F^ + R pX F? (7) 
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and the circular identity 

RpXvn + RpvfiX + RpfiXv = (8) 

the wave equation (6) reduces to the form 

VV^ + R PIU , X F<"> + R\F vp - W v F Xp = (9) 

The Riemann and Ricci curvature coupling terms to the photon polarisation (or spin) 
give rise to the polarisation dependent deviation of photon orbits from null geodesies. 

Photon trajectories are described by the eikonal solutions of the wave equation 
(9) of the form 

F pv = e lS(x) f pv (10) 

where the phase S varies more rapidly in spacetime than the amplitude f pu . The 
wavenumber of the photon trajectories is given by the gradient of the phase 

K p = V P S (11) 

In the geometrical optics approximation where 

V M F a/3 = iK p F a/ 3 (12) 

the wave equation (9) can be written in the form 

- K^K,f uX + u J m - R\f vp + WJ Xp = (13) 

The dispersion relation may be written as 

K 2 = (R p » V J P , - R\U P + R p Jx P ) (14) 
Operating on (14) by V Q , the L.H.S. is 

V a K 2 = 2K^V p K a (15) 
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where we have used the identity V"^ = V ^K 01 which follows from the definition of 

as a gradient. Light rays are defined as the integral curves of the wave vector 
i.e. the curves x^s) for which = K^. Substituting = dX^/ds in (15) we have 
for the L.H.S. 

Using (14), (15) and (16) we obtain the modified geodesic equation for photon tra- 
jectories 

j2va dX@ dX" 1 1 f uX 

IF + r °^^ = 2 V ° ■ f~ - R '^> + ™ " TO (17 > 

To obtain the photon velocity in the curved space, we can use the dispersion 
relation (13) directly. Of the six components of F^ v in equation (13), only three 
are independent owing to the Bianchi identity (5). Choosing the components of the 
electric field vector E,- L = f oi as the independent components we have from the Bianchi 
identity (5) 

KJ :j ■ K;f,„ ■ Ik,/,,, (18) 

Using (18) to substitute for fa in terms of the electric field components fj and fd 
in the wave equation (13) we obtain 

(+K»K, + R\ + R* f ) Sj + (-2X* m + AR' j ol f + R\ 

. \ ( 19 ) 

" lk RJ o k i) foj = 

which is the wave equation obeyed by the three components of the electric field vector 
(in (19) Latin indices are over three spacelike components and repeated indices are 
summed over). To simplify the notation we write (15) as 

{K 2 5{ + 4)f O3 =0 (20) 
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where 

4 = (R o + R e §-) S{ + (-2R°> oi + AR\, + R\ - jfKJc) (21) 

Of the three components of the electric field vector in the electromagnetic waves, 
only the transverse waves are observable in an asymptotically flat space. To obtain 
the equation for the transverse fields we first diagonalise the matrix e\ to give the 
equation for the three normal modes: 

(K 2 + e i )f ot = (i = 1,2,3) (22) 

where e« are the eigenvalues of the e* matrix. The three equations (22) can be sepa- 
rated into the equation for the transverse fields 

ft? = 0} - "S) U (23) 

and the longitudinal fields 

/#> = nS/« (24) 

where n % = are the components of the unit vector along the direction of propaga- 
tion. The wave equation for the transverse photon is 

K 2 (51 - n^n t ) + (5{ - n j m) 

and for the longitudinal photon the wave equation is 

(K 2 + Ci ) (n'm) f oj = (26) 

Expanding Universe (Friedmann- Robert son- Walker Metric): 

For wave propagation in a general homogeneous, isotropic universe described 
by the line element ds 2 = +dt 2 -R 2 {t) (dr 2 + r 2 d6 2 + sin 2 6d<f) 2 ^, the non-zero com- 
ponents of the e* matrix (21) are 

123 „ ( R R2 \ 87rG / 

e!=el = £ S = -2(- + F j=-— (p-3p) (27) 
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fo 3 = 



(25) 



where R is the scale factor, and p and p are the pressure and energy density respec- 
tively. 

The dispersion relation for the transverse photons is 

- kl = ^ (p - 3p) (28) 

and the photon velocity is given by 

(T) duj 1 . . 

= W = 7 W (29) 

Therefore we find that the photons are not superluminal as long as the strong energy 
conditon p > 3p is satisfied. In the radiation dominated era when p = 3p the photon 
velocity is 1. One may turn this argunment around and use the axiom from special 
relativity that photons of any polarisation cannot exceed the flat space velocity c when 
measured in an inertial frame, to show that the strong energy condition p — 3p > 
must be obeyed for any form of matter. Consider a space time region filled with some 
matter with equation of state p = qp. The e*- matrix in that region of spacetime 
is given by e*- = — 8j(8nGp/3) (1 — (3/g)) (assuming the matter distribution to be 
homogenous and isotropic though time dependent). The velocity of photons through 
that medium is given by v = c (1 + (87rGp/3fc?) (1 - q/3)) 1 . Therefore if we assume 
that the photon velocity in an inertial frame does not exceed that flat space velocity 
c, then the equation of state of any form of matter p = qp must obey the strong 
energy condition q > 3. 

Non- Rotating Black Hole Spacetimes (Schwarzschild Metric): 

In the gravitational field of non-rotating blackhole, the non-zero components 
of the 6j matrix (21) are 

i GM 2 n GM 

^ = ±—^1 = 4 = -2-r (so) 



Considering radial trajectories with n = (f, 6, tpj = (1, 0, 0), equation (25) yields the 
wave equations for the transverse fields E 2 ,E 3 : 



K 2 - 2GM/r 3 \/£ 2 

K 2 — 2GM/r 3 { E 3 



(31) 



The dispersion relation yields 



GM 



cu 2 -kf = 2— (32) 



The velocity of propagation of the transverse photon is 

„<<> = = ! < l (33) 

The photon velocity is subluminal and the transverse photons have an effective mass 

/2GM\ 1/2 , x 

m 7 = (— J (34) 

For tangential trajectories, n = (0, 1, 0) and the equations for the transverse fields Ei 
and E 3 yield: 

( K 2 + 4GM/r 3 \ / E 1 \ _ ( . 

{ K 2 — 2GM/r 3 ) { E 3 ) ~ U [6b) 

The dispersion relation is [k 2 + ^j£^J (k 2 — = and the root corresponding 

to the propagating mode is 



2GM 



K 2 =co 2 - k 2 = ^-f- (36) 

and the velocity of wave propagation is given by 

v m = to = ! < ! (37) 

Bh, (l + ^(^))' /2 

In the coordinate frame the dispersion relation obeyed by = (cu, k r kg, k p ) 

obtained from the wave equation (25), for a radial trajectory is given by 

„ 2GM\~ 1 2 / 2GM\ l2 2GM 
1-—) «•-(!-—)*?=— (38) 



Since a; is a constant of motion (the metric being stationary), the wavenumber k r 
changes with r. The redshift of the wavelength A = (g rr ) 1//2 /k r is given by the 
expression 



2(7 .U / 2GM\ ^ 12 
A(r 2 ) = f l-2GM/r 2 \ 1/2 I 1 ~ ~Jyf ^ ~ ~7T/ 



A(n) V 1 " 2 ™/^!/ 1 1 2(7M ^ 2(7M ^ 



(-=r)J 



(39) 



The extra factor in the curly brackets is the correction to the standard redshift formula 
due to the coupling of the Riemann curvature to the wave vector. The correction is 
significant when the photon wavelength is comparable to the Riemann curvature. 

Higher Derivative Gravity: 

Drummond and Hathrell [3] obtained the higher derivative couplings which 
arise from the QED loop corrections to the graviton-photon vertex. The effective 
Lagrangian valid in the frequency range u 2 < a/m 2 e is given by 

C = -iF^ + aRF^Fiv + bR^FTZ 

(40) 

+cR^ aT F^F° T 

with the coefficients a = — JL-^fr = JSt -2 -^ c = — tJtt-^o (where a is the fine 

720 -Km% ' 620 win* ' 720 wm^ \ 

structure constant and m e the electron mass). The equations of motion from (40) are 
given by 

V M F^ + X u = (41) 

where 

X v = a(F» v V u R + RV^F^) 

+2b (F™V U R^ - F^V^R\ - R\V»F°») (42) 

+4c (F^V u R^ aT + R^ aT V u F° T ) = 
Operating on equation (46) by V A and using the Bianchi identity (5) the commutator 
identity (7) and the circular identity (7), we obtain the wave equation in higher 
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derivative gravity: 



+ {RpixvxF^ + R p x F vp — R P U F\ P } 
+ [V V X X - V X X V ] = 



(43) 



In the derivations of [4-12] the terms in the curly bracket do not appear. This is be- 
cause in the derivations, the eikonal approximation is made in the Maxwell's equation 
and therefore the Riemann and Ricci terms which arise on commuting the covariant 
derivatives do not appear. In (43) however we see that in the range of validity of the 
effective action (40), u 2 < a /ml, the Riemann and Ricci terms in the curly brack- 
ets of (48) which are present in Einstein's gravity are larger than the terms in the 
square brackets of (48) which arise from the loop corrections. The photon velocities 
are therefore given by the expression (26) for blackholes, and (39) for Friedmann- 
Robert son- Walker metric, and is subluminal even in the presence of higher derivative 
coupling terms. 
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